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Abstract-  Rigorous  bounds  on  temperature  and  concentration  arc  obtained  lor  a  model  reacting  system 
incorporating  spatial  variation,  reactant  consumption,  and  Arrhenius  kinetics.  Similar  bounds  are  obtained 
tor  the  Semenov  (spatially  uniform)  approximation.  Temperatures  predicted  by  the  Frank  Kamenetskii 
(stationary)  approximation  are  shown  to  be  upper  bounds  on  the  temperature  of  the  complete  time-dependent 
system.  The  latter  supports  numerical  and  approximate  analytic  arguments  that  the  critical  value  of  the 
Frank  Kamenetskii  heat  generation  parameter  in  the  stationary  approximation  is  a  lower  bound  on  analogously 
defined  critical  parameter  values  in  the  complete  system.  The  fundamental  mathematical  tool  is  a  comparison 
result  for  systems  of  parabolic  partial  differential  equations. 


A  comprehensive  model  for  studying  thermal 
explosion  phenomenon  that  incorporates  both 
temporal  and  spatial  variations  in  temperature 
and  reactant  density  is  given  by  the  initial- 
boundary -value  problem 
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(Symbols  are  defined  below  in  Nomenclature. 
Subscripts  denote  partial  derivatives:  for  example. 
Ti *  dTidt*.  The  independent  variable  .v*  may 
be  a  scalar  if  D  is  a  one-dimensional  domain  or  a 
two-  or  three-dimensional  vector  if  0  is  a  plane 
region  or  a  solid.) 

Because  of  the  essential  nonlinearities  and  the 
coupling,  it  is  formidable  to  analyze  the  above 
system  in  its  present  form.  To  make  the  problem 
tractable,  it  is  customary  to  make  various  approxi¬ 
mations.  The  two  commonly  encountered  simpli¬ 
fications  are  the  Semenov  and  Frank  Kamenetskii 
approximations.  In  the  former  case,  spatial 
variation  is  completely  neglected.  The  (hernia! 


dissipation  by  diffusion  through  the  solid  reactant 
is  replaced  by  Newtonian  cooling  at  the  boundary 
of  a  solid  whose  temperature  is  assumed  spatially 
uniform: 

PC i  VT\t*)  ■-  XS(T  T„)  ■  QAVC"  c  *  ltT. 

0)  T„  (5) 

<•'(/*>  AC"  c  "T.  dO)  -  c„.  (6) 

In  the  second  approximation,  one  ignores  reactant 
consumption  while  retaining  thermal  diffusion  in 
a  so-called  stationary  approximation: 
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Many  investigators  have  studied  these  approxi¬ 
mations.  usually  by  imposing  special  conditions  on 
the  parameters  or  by  simplifying  nonlinear  terms: 
for  examples,  see  (Gray  and  Lee.  1967a).  (Mer- 
zhanov  and  A  verson.  1971).  Few  er  hav  e  explored 
relationships  among  these  approximations  and 
the  complete  problem  (I  4):  a  recent  study  that 
emphasizes  connections  between  parameter  groups 
in  the  various  approximations,  especially  when  the 
kinetics  term  v  h:  "r  is  simplified,  is  that  of  Bod- 
dington.  Gray,  and  Wake  (1977). 
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Here  we  shall  take  a  somewhat  different  approach 
to  examining  relationships  among  these  approxi¬ 
mations  and  the  complete  problem.  By  constructing 
solutions  of  appropriate  differential  inequalities, 
we  can.  for  example,  rigorously  establish  that  the 
temperatures  in  (1-4)  or  (5-6)  are  bounded  for  all 
time  or  that  the  temperature  field  defined  by  the 
stationary  approximation  (7-8)  is  larger  at  every 
point  of  the  reacting  solid  than  the  temperature 
field  of  the  complete  problem  (1-4),  regardless  of 
the  time.  We  emphasize  that  these  bounds  are 
mathematically  rigorous  and  that  they  involve  no 
special  restrictions  on  parameters,  approximations 
to  kinetics,  or  the  like. 

To  explain  the  physical  idea  behind  such  bounds, 
suppose  that  the  history  of  reactant  consumption 
at  an  extremely  low  temperature,  say  that  of  the 
surrounding  medium,  is  available.  Since  reactant 
cannot  be  consumed  more  slowly,  this  curve  of 
reactant  density  against  time  must  be  an  upper 
bound  on  the  time  course  of  density  in  the  actual 
reaction.  If  the  reaction  could  now  be  repeated 
w'ith  reactant  being  available  in  amounts  dictated 
by  this  upper  bound  curve,  the  temperature  of  the 
solid  in  this  approximate  reaction  would  be  above 
that  in  the  actual  reaction  because  more  of  the 
exothermic  reactant  would  be  available  at  any 
given  time.  In  other  words,  underestimating 
temperature  would  overestimate  reactant  as  a 
function  of  time.  Overestimating  reactant  would 
then  overestimate  temperature  as  a  function  of 
time. 

After  the  introduction  of  a  set  of  dimensionless 
variables  in  Section  2,  we  derive  the  rigorous 
bounds  on  the  dimensionless  temperature  and 
concentration.  The  relation  of  the  complete 
problem  (1-4)  to  the  approximate  problems  of 
Semenov  and  Frank-Kamenetskii  is  then  con¬ 
sidered  in  Sections  4  and  5,  respectively.  For  the 
sake  of  simplicity  of  exposition,  we  have  considered 
here  a  model  for  a  single  species.  The  techniques 
developed  in  this  paper,  however,  are  capable  of 
handling  multiple  species.  Also,  the  lack  of  a 
diffusion  term  in  Eq.  (2)  (/t*  0)  causes  no 

significant  changes. 


2  DIMENSIONLESS  VARIABLES 

To  nondimensionalize  the  complete  problem  (14) 
and  Frank-Kamenetskii's  stationary  approxi¬ 
mation  (7-8).  we  introduce  the  usual  temperature, 
concentration,  and  length  scales:  see  0.  A.  and  x  in 


Nomenclature.  As  a  time  scale,  we  use  the  charac¬ 
teristic  reaction  time. 

I  it*  =  eK 

the  time  required  to  consume  an  initial  con¬ 
centration  of  c„  in  a  spatially  uniform  reaction 
proceeding  at  the  (constant,  initial)  ambient 
temperature  rate:  />..  l/t/t*  is  the  initial  reaction 
rate.  We  obtain  the  non-dimensionalized  complete 
problem,  denoted  by  (CP): 

a„0,  •  8A"'/(#)  (CP-I) 

X,  —  jiV-A  A >»  f(0)  (CP-2) 

Ht.  ■  Bi  0  0.  A,,  -  0.  xtdD.  t  0.  (CP-3) 

6(x,0)  -  0.  A(.v,  0)  -  I.  x  edD.  (CP-4) 

Proceeding  as  above  but  writing  0E  for  the  scaled 
temperature,  we  obtain  the  non-dimensionalized 
Frank-Kamenetskii  stationary  approximation, 
denoted  by  (FKE): 

0  =  V-i/rE  -  8/(0c)  (FKe-I  ) 

i /it.,,  Biijiz  -  0.  i/itdD.  (FKE-2) 

(If  €  -  0.  then  —  e+  in  (FK„  1).  The 
distinction  between  the  cases  e  0  and  e>0  in  the 
kinetics  term  will  be  important  later.)  The  meaning 
of  i //E  when  (FKE)  exhibits  multiple  solutions  will  be 
clarified  in  Section  5. 

The  parameter  an  is  the  ratio  of  a  characteristic 
diffusion  time  Id*  =  pCtd'lx  to  the  characteristic 
reaction  time  //,•*:  that  is.  ao  is  a  non-dimensional¬ 
ized  characteristic  diffusion  time.  Alternatively. 
aij  is  a  Damkohler  number,  the  ratio  of  reaction 
rate  at  8  =  0  to  the  rate  of  thermal  diffusion.  Of 
course.  8  is  the  usual  Frank-Kamenetskii  para¬ 
meter.  a  non-dimensionalized  heating  rate  or  a 
non-dimensionalized  temperature  increase  when 
diffusion  is  the  only  heat  dissipation  mechanism. 

With  the  same  choice  of  temperature,  reactant, 
and  time  scales,  we  obtain  from  (5-6)  the  non- 
dimensionalized  Semenov  approximation,  denoted 
by  (SA). 

ax0\l) 

A '(/) 


0  >  aX"‘/\0),  0(0)  0.  (SA-I) 

A'"/(0).  A(0)  I.  (SA-2) 
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The  parameters  a, y  and  a  are  analogous  to  "i> 
and  8  in  (CP):  ct.v  =  ix*/iit*  is  a  non-dimensional¬ 
ized  characteristic  time  for  Newtonian  cooling, 
and  a  is  a  non-dimensionalized  temperature 
increase  when  Newtonian  cooling  is  the  only  heat 
dissipation  mechanism.  To  define  /.%•*.  suppose 
T(t)  satisfies  (5)  with  Q  —  0.  If  7\0)  T„  0,  then 
choose  ly*  so  that  T(tx*)  Ta  (T( 0)  T„)/e. 

Note  that  a/{V/xS)  -  8 /(c/2/*-).  Since  the  para¬ 
meter  groups  Vj xS  and  </2/«  are  volume  specific 
heat  loss  terms  for  the  Semenov  approximation 
and  complete  problems,  respectively,  we  see  that 
ihe  heat  generation  term  in  each  problem  has  been 
scaled  against  the  appropriate  dissipation  factor. 

If  the  dissipation  terms  V20  in  (CP-1)  and  6 
in  (SA-I)  are  eliminated,  then  dividing  (CP-2)  into 
(CP-I)  and  (SA-2)  into  (SA-I)  yields,  respectively. 
dBjdX  =  —  hj a u  for  (CP)  and  dBjdX  —  -a/ cry  for 
(SA).  These  quotients  are  equal: 

h/ot)  =  ajay  =  B. 

Since  B=  dd/d\.  which  is  the  temperature 
increase  per  mole  of  reactant  consumed  in  the 
absence  of  heat  loss,  B  is  the  non-dimensionalized 
adiabatic  temperature  rise. 

One  of  the  aims  of  this  paper  is  exploring  the 
effects  of  approximations  of  the  highly  nonlinear 
reaction  terms,  specifically  the  exponential  approxi¬ 
mation.  In  our  notation,  the  exponential  approxi¬ 
mation  corresponds  to  setting  «  =  0  in  the  expres¬ 
sion  for  f(B).  For  example,  in  the  case  of  the 
Frank-Kemenetskii  stationary  approximation. 
(F-K„)  denotes  setting  <  -  0  in  the  kinetics  term 
in  (F-Ke).  thereby  replacing /(ff)  =  exp(0/(l  *-«0)) 
by  e°.  References  to  c  =  0  or  to  infinite  activation 
energy  in  the  following  are  to  the  kinetics  term 
alone,  not  to  any  other  term  or  parameter  appearing 
in  these  nondimensionalized  equations.  Readers 
accustomed  to  asymptotic  analyses  of  these  models 
for  small  c  (large  activation  energy)  should  be 
especially  alert  to  our  treatment  of  «  as  an  in¬ 
dependent  parameter  in  the  kinetics  terms  of  the 
dimensionless  models. 


3  RIGOROUS  BOUNDS 

Here  we  shall  develop  a  technique  for  rigorously 
confirming  the  bounds  on  8(x.t )  and  A(x,  t) 
suggested  by  the  intuitive  arguments  in  the  Intro¬ 
duction.  We  shall  exhibit  systems  of  differential 
inequalities  which  are  bounding  problems  in  the 


sense  that  any  pair  of  functions  satisfying  such  a 
system  of  inequalities  is  a  bound  on  the  pair 
(d(x.  t).  A(.v.  t)). 

To  summarize  the  comparison  theorem  of 
(Chandra  and  Davis.  1979),  which  is  central  to 
these  arguments,  consider  a  given  initial-boundary- 
value  problem  for  a  system  of  parabolic  differential 
equations: 


ut(x.t)  A ,  V  -n  -F{it.  i  ) 

(9) 

i  t(x,  t)  ~  A2V-r  :  G(u,  r) 

(10) 

atu,  ■  h,u  =  0,  CM',,  r  b-,v  —  0,  xedD,  /> 0,  (11) 

u(x.  0)  =  un(x),  r(x,  0)  r„(x).  xtD.  (12) 

We  suppose  that  A,  .  0,  k.,,..0  (to  guarantee  the 
system  is  at  least  weakly  parabolic)  and  that 
at.  bt  0,  ar  '-  br^ 0  if  A/>0  while  at  —  bt  —  0  if 
Af  —  0  (to  guarantee  the  problem  is  neither  over- 
nor  under-determined  by  properly  posed  boundary 
conditions).  We  also  assume  each  smooth  solution 
( u .  c)  depends  continuously  upon  the  parameters 
in  the  differential  equations  and  upon  the  initial 
and  boundary  data.  We  make  no  restrictions  on 
the  monotonicity  of  the  given  nonlinear  terms 
Fand  G. 

However,  we  assume  there  are  functions  F  and  Q 
with  the  properties 

(i)  F(«,  i)^F(u,  c).  G(u,  t')<C(H.  r), 

(ii)  f(u.  r)  is  monotone  nondecreasing  in  r  for 

each  fixed  u  (/>.,  t'<H'  implies  F(u,  tv) 

for  each  u)  and  Q(u,  r)  is  monotone  non¬ 
decreasing  in  it  for  each  fixed  r. 

(Property  (ii)  which  is  essentially  monotonicity  in 
the  off-diagonal  variables  of  f  and  Q.  is  known 
as  quasi-monotonicity  in  the  mathematical  litera¬ 
ture.) 

Finally,  we  suppose  that  (y{x.  t).  r(r,  / ))  is  any 
pair  of  functions  satisfying  the  following  system  of 
differential  inequalities 


IJt(X.  I)' 

A|V*m  •  Fiy.i) 

(13) 

1  r(X.t)< 

kS\  •  Gil<-  f) 

(14) 

•■i.0,  a2i\. 

:  Aji"  0.  xcdD.  t  0, 

(15) 

i-  (a).  | 

[■(.V.  0)'  I„(.v),  xiD. 

(lb) 

156 


.1.  CHANDRA  AND  l\  \V.  DAVIS 


« 


/ 


Comparison  Result:  If  (u.  i)  is  the  solution  of 
(9-12)  and  uj.  r)  is  any  pair  of  functions  satisfying 
(13-16).  then  in.  r)  is  a  lower  hound  on  (u.  r): 

(tt.  i  )■  (u.  i )  (meaning  both  ii(. v.  I )  i/(x.  t )  and 

r(.v. /)■  t(.v.  ; ))  for  all  vc/>  and  all  /  0.)  If /'and 
(i  m  (13)  and  (14)  tire  replaced  by  functions 
satisfying  Properly  (i)  with  all  inequalities  reversed 
ttnd  Property  (ii).  then  solutions  of  <  13  16)  with  all 
inequalities  reversed  are  upper  bounds  on  In.  n. 

We  shall  repeatedly  use  this  Comparison  Result 
by  constructing  appropriate  quasi-monotone  non- 
linearities  derived  from  the  kinetics  terms  of  the 
combustion  problems  in  question  and  by  studying 
various  pairs  of  bounding  functions  satisfying  the 
appropriate  version  of  (13  16). 

Now  consider  the  particular  problem  ICP).  in 
which  u  it.  r  A.  A,  I.  A_.  p.  F  nA'"/( "). 
(i  A'"/( u ).  and  so  on.  To  construct  a  lower 
bounding  problem  whose  nonlinear  terms  satisfy 
Properties  (i)  ttnd  (ii).  we  define 

/  ( u.  A)  ini  | /•( a.  :):  A  r|. 

(J(ti.  A)  inf;(j(r.  A):  t>- 

where  "inf  abbreviates  .inlinum.  or  greatest  lower 
bound.  By  definition.  F  F  and  (J ■  O.  The 
quasi-monotonieity  condition  is  automatically 
satisfied:  F  is  nondeereasing  in  A  (because  the 
inlinum  is  taken  over  smaller  sets  as  A  increases) 
ttnd  (j  is  nondeereasing  in  t>.  We  find  that 
_/<  n.  A)  F(n.  A)  6,\'"f(H)  and  (J(H.  A) 
,\"'f(  j.  )  e1  z\"'. 

1  et  (fj. A)  be  a  pair  of  funetiims  which  satisfy 


I’nL'i ' 

V-6 

.),("'/( 6) 

(17) 

6r  p 

V-A 

<•'  ;  A'" 

(18) 

tj  0. 

A,  0. 

xidl).  t  (>. 

(19) 

0)  o. 

A  (  ,v. 

0)  1.  .Ml). 

(20) 

According  to  the  Comparison  Result.  (17  20)  is  a 
lower  bounding  problem  for  (CP)  in  the  sense  that 
(tj.  A)  (t'(.v.  t ).  A(.v.  / ) )  for  any  pair  (('.  A)  which 
satisfies  (17  20). 

Intuitively,  this  assertion  is  reasonable  because 
(17  20)  models  an  approximate  reaction  which 
starts  with  less,  colder  reactant  (cf.  (20)  and 
(CP-4))  and  proceeds  to  completion  in  a  colder 
environment  into  which  there  is  a  greater  flux  of 
heat  (cf.  (19)  and  (CP-3)).  Reactant  is  consumed 


in  ( IS)  as  if  the  reaction  were  proceeding  at  infinite 
temperature  (/(-/)  e1  :).  Heat  is  released  in 
(17)  by  the  same  mechanism  as  in  the  actual 
reaction,  but  less  reactant  is  available  in  (17)  at 
any  given  time  because  it  is  being  consumed  at  a 
more  rapid,  higher  temperature  rate. 

Note  that  U  0.  a  0  satisfy  (17  20).  thereby 
forcing  0  t'l.v.  / 1.  0  A(.v.  1 1. 

Willi  these  physically  obvious  lower  bounds  in 
hand,  we  can  now  construct  an  upper  bounding 
problem  employing  the  minimum  temperature 

'J  o.  Define 

/ ( ".  a i  sup ;/■(".  ci :  0  a;.  (21) 

0(i'.  ai  sup;(;(.r.  v i :  n  r  (22) 

where  "slip"  abbreviates  supreinum.  or  least  upper 
bound.  We  line)  that  Ii".  i|  /(".  ai  n.\'"/("| 
and  /7(''.  .1)  A"7(0)  A"’, 

l  et  (•(.  >.)  be  a  pair  of  functions  which  satisfy 

V-0  o/. '«/(())  (23) 

7.,  /<V-7.  r"'  (24) 

If  Hi  (1  0.  >.,.  0.  mcI).  t  0.  (25) 

0(.v.0)  0.  7.(.v.  0)  I.  ,v< /).  (26) 

This  system  represents  an  approximate  reaction 
whose  physical  basis  is  the  dual  of  (17  20).  An 
analogous  appeal  to  the  Comparison  Result 
demonstrates  that  (23  26)  is  an  upper  bounding 
problem  for  (CP)  in  the  sense  that  ( "(x.  n. 
A(.v. /))  (0. /.)  for  any  pair  (0.  a)  satisfying  (23  26). 

Similar  arguments  permit  the  construction  of 
bounding  problems  for  (SA).  The  corresponding 
lower  hounding  problem  is 

irxlj'd)  IJ  '<A/'"(6).  6(0)  0.  (27) 

A’)/)  <  1  -A'".  A(  0 )  I.  (28) 

Once  again,  we  may  take  (6.  A)  (0.  0)  to 
conclude  (0.  0)  t "( / 1.  A(/))  for  the  solution 
( 0(1 ).  A(mi  of  (SA).  (  sing  (21  22).  we  then  define 
an  upper  bounding  problem  for  (SA): 

ov'V(M  0  .i7."'/(0).  0(0)  0.  (29) 

7.0)  7.'".  7.(0)  I.  (30) 
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To  exhibit  a  bound  on  the  solution  of  (CP)  or 
(SA),  we  need  only  find  a  pair  of  functions  satis¬ 
fying  the  appropriate  bounding  problem.  Since 
these  bounding  problems  are  defined  by  differential 
inequalities,  one  may.  in  principle,  construct  a 
number  of  different  bounding  pairs  from  a  single 
bounding  problem. 

4  BOUNDS  FOR  THE  SEMENOV 

APPROXIMATION  AND  THE  COMPLETE 
PROBLEM 

Without  making  the  assumption  «  0.  we  shall 

obtain  spatially  uniform  bounds  on  temperature 
and  reactant  concentration  for  both  (SA)  and 
(CP)  by  studying  a  pair  of  ordinary  differential 
equations  whose  solution  satisfies  an  upper 
bounding  problem  for  (SA)  and  (CP).  Roughly, 
these  bounds  show  that  bounded  kinetics  (/.<’.. 
finite  activation  energy  or  «  0)  imply  bounded 
temperatures.  When  <  0.  in  which  case  the 

kinetics  term  is  unbounded  due  to  infinite  acti¬ 
vation  energy  ,  these  techniques  provide  a  bound  on 
induction  time  based  on  reactant  consumption 
occurring  at  ambient  temperature  (H  0). 


Let  (0(i).  /.(/))  denote  the  solution  of 

0(; )  fl/.'"(/)/(0(/)).  0(0)  0 

(31) 

'/.'(/)  /.'"(M.  /.((»)  1. 

(32) 

If  B  denotes  <r/,.  then  (0(/>.  >.(/))  evidently 
satisfies  (23  26).  an  upper  bounding  problem  for 
(CP).  If  B  denotes  «  «\.  then  (<)(;).  >.(;))  satisfies 
(29  30).  a  corresponding  upper  bounding  problem 
for  (SA).  That  is.  (0(m.  /.(/>(  (i>.  A)  whether 
(ii.  A)  satisfies  (SA)  or  (CP). 

To  obtain  an  expression  lor  0(/).  substitute  (32) 
in  (31).  separate  variables,  and  integrate: 

"in 

1  Jn  /un  Hi  1  uni. 

n 

(33) 

Solving  (32)  independently  yields 

/.(/)  <  '.  in  1. 

(34a) 

/.(/ )  [1  (m  1)/)'  1  in  ■■  1 

(34b) 

(When  0  m  I  in  (34b).  />(/)  takes  the  form  given 
for  ()•  r  III  »i)and/.(r)  0  for  1(1  m)  /.) 


Since  lint  /.(/ )  0,  lim  U(t )  ~  ()( x )  must  satisfy 

(  -•*  (~*y. 

*k  *  ) 

|  duif(u)  =  B.  (35) 

0* 

Since  _()'(/)  0  in  (31).  )(/)  is  increasing  and 
0 (/)•  0(x)  for  all  l.  Hence.  0(x)  is  an  upper 
bouitd  on  the  maximum  temperature  in  either 
(SA)  or  (CP). 

From  (35).  we  sec  that  this  bound  0(  x )  on  the 
maximum  temperature  depends  only  the  dimen¬ 
sionless  adiabatic  temperature  rise  B.  Neither 
dissipation  nor  Newtonian  cooling  is  included. 
The  problem  (31  32)  defining  (fi.  /.)  is  similarly 
restricted  as  well;  (31)  could  be  thought  of  as  an 
approximation  to  (SA-1)  or  (CP-I)  in  which 
dissipation  has  been  eliminated  by  discarding  the 
term  l>  or  Since  the  spatially  uniform 

function  0(/)  satisfies  the  boundary  condition 
(CP-3)  with  Bi  0.  <)(/)  may  be  thought  of  as  the 
temperature  of  a  completely  adiabatic  reaction. 

The  supply  of  reactant  to  this  adiabatic  reaction 
is  governed  by  (32).  which  is  the  rate  equation  at 
ii  0.  That  is.  this  hypothetical  adiabatic  reaction 
has  more  reactant  available  at  any  time  than  does 
(SA;  or  (CP)  because  it  loses  reactant  at  the 
slowest,  ambient  temperature  rate. 

These  remarks  provide  intuitive  confirmation 
of  the  rigorous  result  that  temperature  and  reactant 
concentration  in  either  (SA)  or  (CP)  are  bounded 
above  by  (0(r ).  /.(,'))  and  that  the  maximum 
temperature  in  either  case  is  bounded  by  0(x)). 
Further,  they  warn  us  that  fi(x)  may  substantially 
overestimate  the  true  maximum  temperature. 

To  show  that  0(x)  is  always  finite,  note  that 
fin)  /(  /  )  t’Vl.  Then  (35)  yields 

A  ' ) 

B  |  c  1  ’  (In  =  e  1  ’  fi(  x ), 
n 

or 

<1(  xKfle'e.  (36) 

Boddington.  Gray,  and  Wake  (1977;  Eq.  (14)) 
have  obtained  a  similar  upper  bound.  B  n\c'/s  in 
our  notation,  in  the  uniform  temperature  case 
without  reactant  consumption  but  with  dissipation 
(».<-..  (SA-I)  alone  with  A  s  I).  We  now  see  that  a 
bound  of  this  type  applies  as  well  to  both  the 
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uniform  temperature  case  with  reactant  con¬ 
sumption  (SA).  as  one  might  expect  on  physical 
grounds,  and  to  the  distributed  temperature  case 
with  reactant  consumption  (CP),  where  the  physical 
relation  to  the  uniform  temperature  case  is  perhaps 
less  obvious. 

Put  simply,  (36)  implies  that  the  temperature  in 
both  (SA)  and  (C'P)  remains  finite  for  all  time  so 
long  as  the  activation  energy  is  finite  («  0). 

Unfortunately,  (35)  cannot  provide  tighter 
bounds  on  the  maximum  temperature.  Since 
/(«).  I.  (35)  yields 


ft  <  ) 


In  physical  terms,  a  bounding  problem  that 
neglects  all  dissipation  will  exhibit  maximum 
temperatures  at  least  as  great  as  the  adiabatic 
temperature  rise. 

We  may.  however,  use  the  upper  bound  (36) 
to  define  an  improved  lower  bounding  problem  for 
either  (SA)  or  (CP)  by  replacing  /(  x )  c' •  in 
(18)  or  (28)  by  /{Be'1).  If  A(/ )  now  satisfies 

A</ )  A  "‘(t)f(Bc'z),  ^(0)  I. 

then  the  pair  (0.  A(/ ))  satisfies  either  of  these 
improved  lower  bounding  problems.  We  calculate 

A( t )  exp[  J\Bc'i')t],  m  =■■=--  I. 

A(0  (I  (/;/  l).f(Bc'  *)/]’■■  (i  «»,  m  /  I. 

(If  0  w  I  in  the  last  expression,  then  A(/)  takes 
the  form  given  for  0  r  [(1  ;»)/(fic’ c)]  1  and 

A(0  0  for  larger  values  of  /.) 

Hence,  if  Ml)  is  defined  by  (34)  and  (II.  A)  is  a 
solution  of  either  (SA)  or  (CP), 

(O.A(O)  in.  A)-  (f)(/),  7.(0)  (Bc'^MD).  (37) 

Since  the  dimensionless  time  t  has  been  scaled 
against  In*,  the  characteristic  reaction  time,  (37) 
suggests  that,  in  dimensional  time,  the  time 
constant  for  reactant  consumption  ranges  between 
I J  (Be1  z)t  n*  and  I  in*. 

If  we  make  the  exponential  approximation 
j\0)  *  c".  then  this  unbounded  kinetics  term 
permits  0(/ )  in  the  adiabatic  bounding  problem 
(31  32)  to  grow  without  bound.  To  calculate  the 


induction  lime  at  which  0  -  x,  write  c"  in 
place  of /(</)  in  (33)  and  integrate: 

I  c  lhn  B(  1  7.(0). 

If  <)(/',)  x .  then  I  fi(l  >.(/,)).  and 
/,  /.‘(I  I  fi). 

where  7.  1  denotes  the  inverse  of  the  function  /. 
defined  by  (34).  Specifically,  we  have 

f.  Inf  I  - 1  fi).  m  I. 

(l-l/fl)1  '"-I 

ti  =  — — - hi  yt  1 . 

in  -  I 

(Note  that  fi  I :  extremes  are  4  fi  1.000  while  a 
typical  value  for  spontaneous  ignition  of  wood  is 
fl  ^  30  (Tyler  and  Wesley.  1967).  (Gray  and  Lee, 
1967b).)  In  either  case,  an  easy  argument  shows 
fi  I  S  or.  in  dimensional  variables. 

fl*  0d;4:  hr,  B  f.„ i*. 

where  o.i*  is  the  induction  time  for  the  spatially 
uniform  adiabatic  reaction  with  infinite  activation 
energy  and  no  reactant  consumption.  Of  course, 
an  induction  time  such  as  fr  computed  with 
reactant  consumption  should  be  larger  than 
an  induction  time  computed  without  reactant 
consumption,  as  illustrated. 

If  0(0  continues  to  denote  the  solution  of  (31) 
with  e  0.  then  our  same  bounding  arguments 
can  show  that  0(0  is  an  upper  bound  on  0  in  either 
(SA)  or  (CP)  with  c  0  (or  e  0  as  well).  Hence, 
neither  of  these  problems  can  exhibit  infinite 
temperature  before  0(0  x .  That  is.  the  in¬ 
duction  period  U  for  cither  (SA)  or  (CP)  at  infinite 
activation  energy  must  satisfy 

U  .  '  U  .  [  B  -  tail. 

Hence,  fi  is  a  lower  bound  on  the  induction  time 
in  either  (SA)  or  (CP)  with  c  0. 


5  THL  RELATION  BETWEEN  THE 

STATIONARY  APPROXIMATION  AND 
THE  COMPLETE  PROBLEM 

The  basic  result  of  this  section  is  the  following:  if 
Mi)  is  defined  by  (34),  then 


/ 
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(0u(.v).  /.(')>  (0S(.V). /.(/)).  (6(.v.  ( ).  A(.v.  I )),  (38) 


Arif  Am*, 


where  0£  is  the  stationary  temperature  of  (I  K.-l 
and  0  is  the  exact  temperature  of  (CP).  If  A  in 
(FK„)  is  supraeritical.  then  >/>ii(.y )  does  not  exist  and 
only  the  second  inequality  is  valid. 

In  other  words,  either  of  the  stationary  tempera¬ 
tures  i/i„(.v)  or  i/<j(.v)  hounds  the  true  temperature 
(l(.v.  I )  given  by  (CP). 

The  first  inequality  in  (38)  is  a  consequence  of 
j\0)  and  the  scalar  maximum  principle  (Prottcr 
and  Weinberger.  1967.  p.64).  To  repeat  the 
argument  in  slightly  different  terms,  the  relation 
e'F  •/((()  makes  (FK„)  tin  upper  bounding  problem 
for  (FK£). 

To  establish  the  second  inequality  in  (38).  we 
first  note  that  direct  calculation  in  (34)  shows 
>.(/)•  I  while  the  maximum  principle  for  elliptic 
boundary  value  problems  (Plotter  and  Weinberger. 
1967.  p.64)  forces  (M  v)  0.  The  former  y  ields 

<r/>'/<c.|(.V)  0  V-’i/zj  iV/(i/j:)  V-i/i.  <3>,'"/(i/<£) 

(39) 


a  phenomena  observed  in  numerical  computations 
by  Tyler  and  Wesley  ( 1967)  and  Parks  ( 1961 )  and  in 
various  approximate  analytic  studies  by  Frank- 
Kamenetskii  (1945)  (corrected  in  (Gray  and  Lee, 
I9t»7b)).  Gray  and  Lee  (1967b),  Thomas  (1961), 
Boddington,  Gray,  and  Wake  (1977).  c ft  al.  For  a 
specific  comparison  in  the  spatially  uniform  case, 
see  (Kassoy  and  Lilian.  1978.  Lq.  (5.15)).  That  is. 
the  critical  value  of  the  Frank-Kamenetskii 
parameter  appears  to  be  a  lower  bound  on  critical 
values  for  (FK£)  or  (CP)  defined  by  sharp  increases 
in  temperature. 

Just  as  in  the  previous  section,  the  upper  bounds 
on  temperature  of  (38)  lead  to  new  lower  bounds  on 
reactant  concentration.  However,  these  bounds 
exhibit  spatial  variation  as  well. 

Define  VJ(\.  i ).  A(.v.  i ))  by 


,  V-6 

(40) 

/‘T-A  — A"'/(<//£(.v)). 

(41) 

while  the  latter  ensures  that  i bz  satisfies  the  initial 
condition  (26)  in  the  upper  bounding  problem. 

Comparing  (39)  and  the  delining  equation  (32) 
for  /.(/)  with  (23)  and  (24).  respectively,  shows  t hat 
( 0 £(.v),  /.(/))  is  a  solution  of  the  upper  bounding 
problem  (23-26)  for  (CP).  Hence.  ('C(.v). 
7.(0)  (O(.v.  t),  A(.v.  /)). 

If  b  is  chosen  so  thill  (I  K£).  e  0.  has  several 
solutions  of  which  only  one  is  stable,  then  i/»£(.v)  is 
the  smallest  non-negative  stable  solution  satisfying 
(25).  for  stability  is  the  manifestation  of  the  con¬ 
tinuous  dependence  condition  of  the  Comparison 
Result  (Chandra  and  Davis,  1979).  upon  which 
rests  the  validity  of  our  bounding  claims.  If  a  is 
in  a  range  in  which  (FKt)  has  several  stable 
solutions,  then  i//£(.v)  may  denote  any  one  of  these, 
including,  of  course,  the  smallest. 

Suppose  bt  and  An-  are  “critical”  values  of  a  in 
(FK£)  and  (CP),  respectively,  defined  by  some 
reasonable  criteria  involving,  say.  a  sharply 
increased  maximum  temperature.  (Recall  that 
results  of  the  previous  section  guarantee  that  the 
temperature  in  (FK£).  <  0.  or  (CP)  is  never 
infinite.)  Let  f>,rit  be  the  usual  Frank  Kamenet- 
skii  critical  value  for  (FK„)  “beyond  which  this 
stationary  distribution  ceases  to  be  possible” 
(Frank  Kamcnetskii,  1969,  p.345).  Then  the 
bounds  of  (38)  suggest 


())  0.  A(.v.O)  I.  xel)  (42) 

ij,  Bit!  0.  A,,  -  0  \(dl).  i  •(>.  (43) 

Then  (L1-  A)  evidently  satisfies  an  improved  lower 
bounding  problem  obtained  from  (17-20)  by 
replacing/)-/)  :  in  ()8)  by  /(i/;£(.v)).  Hence. 
U'(.v.  I).  A(.v.  ())  ((('(.v.  I).  A(.v.  /)).  where  (0,  A) 
stilislics  (40  43). 

For  a  given  reactant  geometry,  one  could,  in 
principle,  compute  </<t(.v)  and  solve  (41-42)  for 
A(.v.  / ).  Then  the  expression  for  A(.v.  I)  could  be 
used  in  (40)  to  obtain  a  single  differential  equation 
for  a  non-trivial.  spatially  varying  lower  bound  on 
temperature.  The  expression  for  A(.v.  I)  would  also 
provide  a  spatially  varying  lower  bound  on  the 
time  constant  for  reactant  consumption:  see  the 
remarks  below  (37). 


6  CONCLUSIONS 


Bounds  have  been  constructed  on  temperature  and 
reactant  concentration  in  the  Semenov  approxi¬ 
mation  (SA)  (spatially  homogeneous  with  reactant 
consumption)  and  in  the  complete  problem  (CP) 


(spatially  distributed  with  reactant  consumption) 
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|l>0 

These  bounds  show  (hat  temperature  and  con¬ 
centration  are  never  negative  and  (hat  the  (dimen¬ 
sional)  temperature  in  either  case  is  hounded  above 

by 

EQcae*  "'C,  ,>CvRT,r  (44) 

so  long  as  activation  energy  is  finite.  The  (dimen¬ 
sional)  time  constant  Tor  reactant  consumption  in 
either  case  ranges  between 

,4c,,'"  ’  ltr«f\c„QEcli  ,>Cv  RTr) 
and  ,4r„"'  1  <■'•'  ltT». 

The  temperature  bound  (44)  exceeds  the  adia¬ 
batic  temperature  cuQbC which  is  the  maximum 
physically  allowable  temperature  rise,  because  the 
approximate  model  (31-32)  from  which  it  was 
derived  overestimates  fuel  supply.  This  is  one 
example  of  a  common  occurrence  with  the  bounding 
techniques  used  here:  one  can  often  identify  the 
physical  basis  for  a  mathematically  rigorous  bound. 
Such  a  physical  observation  can  then  suggest  a 
mathematical  approach  to  improved  bounds. 
For  another  example,  see  (40-43)  and  the  next  ti¬ 
the  last  paragraph  of  the  Introduction. 

In  the  exponential  approximation  (exp[(<  ( I  «»)] 
'  c").  which  corresponds  to  infinite  activation 
energy  in  the  kinetics  term,  temperature  predicted 
by  the  models  is  not  necessarily  bounded.  In  this 
case,  the  (dimensional)  induction  period  is  no 
shorter  than 

-  Ac„'"  '  <’  "r«  ln(  I  ,,CvRT„-  h:Qc„). 

m  =  I . 

.4c,,’"'1  e  /;  -  I )  1 

[(I  / iCcRTa-'EQc „)’  I),  in  •  I. 

Recalling  our  convention  that  references  to  infinite 
activation  energy  («  0)  arc  to  the  kinetics  term 

alone,  we  note  that  infinite  dimensionless  tempera¬ 
ture  ft  corresponds  to  infinite  dimensional  tempera¬ 
ture  T  and  conversely. 

The  temperature  in  the  complete  problem  (C  P), 
the  spatially  distributed  reaction  with  reactant 
consumption,  is  always  bounded  above  by  the 
temperature  field  of  the  stationary  approximation 
with  finite  activation  energy,  (FK£).  This  stationary 
temperature  at  finite  activation  energy  is  in  turn 
bounded  above  by  the  temperature  in  Frank  - 
Kcmenetskii's  stationary  approximation  with 
infinite  activation  energy.  (FK„).  Consequently, 
one  anticipates  that  critical  behavior  defined  by 
sharp  increases  in  temperature  must  occur  in  the 


opposite  order:  the  highest  temperature  reaction 
(FK„)  can  exhibit  critical  behavior  while  the  other 
reactions  (FK£)  and  (CP)  can  still  be  stibcritical 
at  the  same  parameter  values,  and  so  on. 

In  particular,  the  critical  parameter  value 
<Vr,i  of  Frank- Kcmenetskii's  infinite  activation 
energy,  stationary  approximation  appears  to  be  a 
lower  bound  on  analogously  defined  critical 
parameter  values  for  both  the  finite  activation 
energy,  stationary  approximation  (FK£)  and  the 
complete,  spatially  distributed  problem  with 
reactant  consumption  (CP).  This  observation  is 
based  on  the  rigorous  ordering  derived  above  of 
the  temperature  fields  in  the  three  problems 
(FK„).  (FK£).  and  (CP),  thereby  providing  an 
alternate  view  of  the  mechanism  which  orders 
critical  parameters  in  the  three  cases.  The  three 
problems  do  not  appear  to  have  been  considered 
simultaneously  in  this  light  by  other  investigators. 

The  results  presented  here  do  not  depend  on 
special  restrictions  on  parameters  or  on  approxi¬ 
mations  to  nonlinear  terms. 


NOMF.NT  LATTJRF 

.4  pre-exponential  factor  in  Arrhenius  kinetics 

term  (mofm8)1  ms) 

B  EQc„lpC\  RT<r.  dimensionless  adiabatic 

temperature  rise 
Bi  x‘t  *■  Biot  number 

<  reactant  concentration  (mol/m3) 

c„  initial  reactant  concentration  (mol/m3) 

Cr  reactant  specific  heat  capacity  at  constant 

volume  (J/K  kg) 

</  characteristic  length  of  region  enclosing 

reactant  ( m ) 

I)  region  enclosing  reactant 

<  /)  boundary  of  region  enclosing  reactant 

/.  activation  energy  (J/mol) 

/{»)  exp(((  (l  :  (H)).  dimensionless  kinetics 
factor 

w  reaction  order  (m  0) 

Q  exolhermicity  per  mole  (J,  mol ) 

R  universal  gas  constant  (J/K  mol) 

S  surface  area  of  region  enclosing  reactant 

(in3) 

i*  time  (s) 

i  I*  in  4 c„"'  1  I*  ch:  ,ST-.  dimensionless 

time 

/„„*  in*  B  nCcRTa-'c*1^-  F.(>4 <s) 

induction  time  for  spatially  uniform 
adiabatic  reaction 
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